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We present numerical A'^-body simulation studies of large-scale structure formation. The main 
purpose of these studies is to analyze the several models of dark matter and the role they played in 
the process of large-scale structure formation. We analyze the standard and more successful case, 
i.e., the cold dark matter with cosmological constant (ACDM). We compare the results of this model 
with the corresponding results of other alternative models, in particular, the models that can be built 
from the Newtonian limit of alternative theories of gravity like scalar-tensor theories. An specific 
model is the one that considers that the scalar field is non-minimally coupled to the Ricci scalar in 
the Einstein-Hilbert Lagrangian that gives, in the Newtonian limit an effective gravitational force 
that is given by two contributions: the standard Newtonian potential plus a Yukawa potential that 
comes from a massive scalar field. Comparisons of the models are done by analyzing the snapshots 
of the A''-body system at z=0 for several values of the SF parameters. 

PACS numbers; 95.30.Sf; 95.35. +d; 98.65.-r; 98.65.Dx 



I. INTRODUCTION 



In this work we present some results about the role scalar fields (SF) play in cosmological simulations, in particular 
on the process of large scale structure formation. The main goal of this work is to study the large scale structure 
formation where the usual approach is that the evolution of the initial primordial fluctuation energy density fields 
evolve following Newtonian mechanics in an expanding background The force between particles are the standard 
Newtonian gravitational force. We will see that we can introduce SF by adding a term in this force. This force will 
be of Yukawa type with two parameters (a, A)j3|. We have been studying, in the past years, the effects of this kind of 
force on some astrophysical phenomena [2, y, |j, |5|. The Yukawa force comes as a Newtonian limit of a scalar-tensor 
theory (STT) with the SF non-minimally coupled to gravitation Q. It is our purpose to find the role these scalar fields 
play on the large scale structure formation processes. We start by discussing the standard ACDM model and the 
general approach in A^-body simulations (See Bertschinger^] for details). Then, we present the modifications we need 
to do to consider the effects of a static SF and we show the results of this theory for the cosmological concordance 
model of a ACDM universe 0. To perfor m the simulations we have modified a standard serial treecode the author 
has developed 8] and the Gadget 1 Q (see also http://www.astro.inin.mx/mar) in order to take into account the 
contribution of the Yukawa potential. 



II. EVOLUTION EQUATIONS FOR A ACMD UNIVERSE 
A. General Scalar-tensor theory and its Newtonian limit 

Let us consider the Einstein field equations of a typical STTfloj 
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for the metric 5^1^ and for the massive SF (p we have 
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where ()' = Here T^u is the energy-momentum tensor with trace T, and V{(f>) are in general arbitrary 

hmctions that gobern kinetic and potential contribution of the SF. The potential contribution, V{(t>), provides mass 
to the SF, denoted here by msF- 
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The study of large-scale formation in the universe is greatly simplified by the fact that a limiting approximation 
of general relativity, Newtonian mechanics, applies in a region small compared to the Hubble length cH^^ {cH^^^ « 
3000/i~^ Mpc, where c is the speed of light, Hq = lOO/i km/s/Mpc, is Hubble's constant and h w (0.5 — 1)), and 
large compared to the Schwarzschild radii of any collapsed objects. The rest of the universe affect the region only 
through a tidal field. The length scale cHq^ is of the order of the largest scales currently accessible in cosmological 
observations and H^^ « 10^°/i~^ yr characterizes the evolutionary time scale of the universe. 

Therefore, in the present study, we need to consider the influence of SF in the limit of a static STT, and then we 
need to describe the theory in its Newtonian approximation, that is, where gravity and the SF are weak (and time 
independent) and velocities of dark matter particles are non-relativistic. We expect to have small deviations of the 
SF around the background field, defined here as {(j)) and can be understood as the scalar field beyond all matter. If 
one defines the perturbations (f> ={(/))+ (p and (^^j^ = Ty^^ + h^^, where ry^i^ is the Minkowski metric, the Newtonian 
approximation gives [3] 

Roo - Jv^/ioo - ^4^P - l^^4> , (3) 
2 1 + a 2 

V^^ - mlpct) = -%TTap , (4) 

we have set (0) = (1 + a)/GN and a = 1/(3 + 2tj). In the above expansion we have set the cosmological constant 
term equal to zero, since on galactic scales its influence should be negligible. We only consider the influence of dark 
matter due to the boson fleld of mass msp governed by Eq. (|4]), that is the modified Helmholtz equation. However, 
at cosmological scales we do take into account the cosmological constant contribution, see below. Equations ^ and 
(jl]) represent the Newtonian limit of the STT with arbitrary potential V{(j>) and function u}{4>) that where Taylor 
expanded around (4>). The resulting equations are then distinguished by the constants Gn, ct, and A — hp/mspc. 
Here hp is Planck's constant. 

The next step is to find solutions for this new Newtonian potential given a density profile, that is, to find the 
so-called potential-density pairs. General solutions to Eqs. ([3]) and ([4]) can be found in terms of the corresponding 
Green functions, and the new Newtonian potential is0, 0] 

2 l + aj |r-r^ 

-„i^/,. /('f-"7''\ B.C. (5) 

The first term of Eq. ([5|), is the contribution of the usual Newtonian gravitation (without SF), while information 
about the SF is contained in the second term, that is, arising from the influence function determined by the modified 
Helmholtz Green function, where the coupling uj (a) enters as part of a source factor. 



B. Cosmological evolution equations using a static STT 



To simulate cosmological systems, the expansion of the universe has to be taken into account. Also, to determine 
the nature of the cosmological model we need to determine the composition of the universe, i. e., we need to give the 
values of flj for each component i, taking into account in this way all forms of energy densities that exist at present. 
If a particular kind of energy density is described by an equation of state of the form p = wp, where p is the pressure 
and w is a constant, then the equation for energy conservation in an expanding background, d{pa^) = —pdla^), can 
be integrated to give p oc a"^^^^™-*. Then, the Friedmann equation for the expansion factor a{t) is written as 
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where Wi characterizes equation of state of specie i. The most familiar forms of energy densities are those due to 
pressureless matter with = (that is, nonrelativistic matter with rest-mass-energy density pc^ dominating over 
the kinetic-energy density pw^/2) and radiation with Wi — 1/3. The density parameter contributed today by visible, 
nonrelativistic, baryonic matter in the universe is k, (0.01 — 0.2) and the density parameter that is due to radiation 
is ^p « 2 X 10~^. In this work we will consider a model with only two energy density contribution. One which is a 
pressureless and nonbaryonic dark matter with ^dm ~ 0.3 that does not couple with radiation. Other, that will be a 
cosmological constant contribution f^A ~ 0.7 with and equation of state p = —p. The above equation for a[t) becomes 
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Here, we employ a cosmological model with a static SF which is consistent with the Newtonian Umit given by Eq. 
([5]). Thus, the scale factor, a{t), is given by the foUowing Friedman model. 



1 + a V 1 + a 



(8) 



where H — a/a, flmo and Oaq are the matter and energy density evaluated at present, respectively. We notice that the 
source of the cosmic evolution is deviated by the term 1 + a when compared to the standard Friedman-Lemaitre model. 
Therefore, it is convenient to define a new density parameter by Sll"-* = + a). This new density parameter is 

such that flm^ + fi^"^ ~ which implies a flat universe, and this shall be assumed in our following computations, 
where we consider {flm\^^j^'') = (0.3,0.7). For positive values of a, a flat cosmological model demands to have a 
factor (1 + 0;) more energy content (fim and Q\) than in standard cosmology. On the other hand, for negative values 
of a one needs a factor (1 + a) less ft,n and JIa to have a flat universe. To be consistent with the CMB spectrum and 
structure formation numerical experiments, cosmological constraints must be applied on a in order for it to be within 
the range (-1,1) [H [H d [ll . 

C. The A'^-Body problem for dark matter 

The Vlasov-Poisson equation in an expanding universe describes the evolution of the six-dimensional, one-particle 
distribution function, /(x, p). The Vlasov equation is, 

l + ^-l^— V$.(x).|^.0 (9) 
ot ma'' ox op 

where x is the comoving coordinate, p = ma^x, m is the particle mass, and <I>Ar is the self-consistent gravitational 
potential given by the Poisson equation, 

V2$w(x) = 4^Gjv a' [p(x) - pf,(i)]] (10) 



where pi, is the background mass density. Eqs. ^ and (|T0|) form the Vlasov-Poisson equation, constitutes a collision- 
less, mean-field approximation to the evolution of the full A^-body distribution. An A^-body code attempts to solve 
Eqs. ([5]) and pop by representing the one-particle distribution function as 



N 

/(x,p) = ^<5(x-x05(p-p.) (11) 

1=1 

Substitution of pT|) in the Vlasov-Poisson system of equations yields the exact Newton's equations for a system of N 
gravitating particles. See Ref. [1] for details. 

In the Newtonian limit of STT of gravity, the Newtonian motion equation for a particle i is written as[20j 

where the sum includes all periodic images of particle j, and Fspif, ct, A) is 

FsF(r,a,A) = 1 + a (l + ^) e-'-/^ (13) 

which, for small distances compared to A, is Fspii" < A, a. A) « 1 + a (l + and, for long distances, is Fspir > 
A, a, A) « 1, as in Newtonian physics. 

III. RESULTS 

In this section, we present results of cosmological simulations of a ACDM universe with and without SF contribution. 
We use 256'^ particles in box 256 M pc size. We have studied in the past a ACDM model in a smaller box and 
with less resolution than the present case[20|, the ACDM case that comes with Gadget 1.0[19], and the Santa Barbara 



cluster [21 1 
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FIG. 1: x-y snapshots at z = of a ACDM universe without and with SF. See text for details. 



The initial linear power spectrum was generated using the fitting formula by Klypin & Holtzman[l5| for the 
transfer function. This formula is a slight variation of the common BBKS fit[16|. It includes effects from baryon 
suppression but no baryonic oscillations. We use the standard Zel'dovich approximation [l^ to provide the initial 
256'^ particles displacement off a uniform grid and to assign their initial velocities in a 256h~^ Mpc box. The 
starting redshift is = 50 and we choose the following cosmology: floM — 0.314 (where ftoM includes cold dark 
matter and baryons), ^Ib = 0.044, SIa = 0.686, Hq = 71 km/s/Mpc, as = 0.84, and n = 0.99. Particle masses 
are in the order of 1.0 x 10^'^ Mq. The individual softening length was 20 kpc/h. This choice of softening length 
is consistent with the mass resolution set by the number of particles. All these values are in concordance with 
measurements of cosmological parameters by WMAPpJj]. The initial condition is in the^Cosmic Data Bank web page 
(http://t8web.lcaii.gov/people/heitmann/test3.htinl). See Heitmann et al. 2005[7j for more details. 

We now present the results for the ACDM model previously described. Because the visible component is the smaller 
one and given our interest to test the consequences of including a SF contribution to the evolution equations, our 
model excludes gas particles, but all its mass has been added to the dark matter. We restrict the values of a to the 
interval (—1, 1) [111 [H, [Hi 13 and use A = 1, 5, 10, 20 Mpc/h, since these values sweep the scale lengths present in 
the simulations. In Fig. 1 we show x-y snapshots at redshift z — oi our ACDM model. Fig. 1 (a) presents the 
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standard case without SF, i.e., the interaction between bodies is through the standard Newtonian potential. In (b), 
(c), (d), and (e) we show the case with a = 1, and A = 1,5, 10,20 Mpc/ft., respectively. In (f) and (g) A = 5 Mpc/h 
and a = —1/2 and —1/4, respectively. One notes clearly how the SF modifies the matter structure of the system. 
The most dramatic cases are (e) and (f) where we have used a — 1 and A = 20 Mpc//i, and a — —1/2 and A = 5 
Mpc/h, respectively. 

We now analyze the general effect that the constant a has on the dynamics. The role of a in our approach is as 
follows. On one hand, to construct a flat model we have set the condition flm"^ + fJ^'' = 1, which implies having 
(1 + a) times the energy content of the standard ACDM model. This essentially means that we have an increment by 
a factor of (1 + a) times the amount of matter, for positive values of a, or a reduction of the same factor for negative 
values of a. Increasing or reducing this amount of matter affects the matter term on the r.h.s. of the equation of 
motion p^ . but the amount affected cancels out with the term (l + a) in the denominator of (|12p stemming from the 
new Newtonian potential. On the other hand, the factor Fsf augments (diminishes) for positive (negative) values of 
a for small distances compared to A, resulting in more (less) structure formation for positive (negative) values of a 
compared to the ACDM model. For r ^ A the dynamics is essentially Newtonian. 

Therefore, for cases in which we use A = 5 Mpc/h we have the following. In the case of (c), for r <^ X, the effective 
gravitational pull has been augmented by a factor of 2, in contrast to case (f) where it has diminished by a factor 
of 1/2; in model (g) the pull diminishes only by a factor of 3/4. That is why one observes for r < A more structure 
formation in (c), less in (f), and lesser in model (g). The effect is then, for a growing positive a, to speed up the 
growth of perturbations, then of halos and then of clusters, whereas negative a values (a — > —1) tend to slow down the 
growth. Whereas in models (b), (c), (d) and (e) where we keep a = 1 and take A = 1, 5, 10, 20 Mpc/h, we observe less 
structure in case (b) to more structure in case (e), passing for intermediate structure formation cases. In spite of that 
the effective gravitational constant has been augmented by a factor of 2 the importance of the Yukawa contribution 
is very small for distances r 3> A. That is way we observe this behavior. 

IV. CONCLUSIONS 

The theoretical scheme we have used is compatible with local observations because we have defined the background 
field constant < (f> >— G]^^{1 + a). A direct consequence of the approach is that the amount of matter (energy) has 
to be increased for positive values of a and diminished for negative values of a with respect to the standard ACDM 
model in order to have a flat cosmological model. Quantitatively, our model demands to have Jl/(1 + a) = 1 and 
this changes the amount of dark matter and energy of the model for a flat cosmological model, as assumed. The 
general gravitational effect is that the interaction including the SF changes by a factor Fsrir, a. A) ~ 1 + a (l + 
for r < A in comparison with the Newtonian case. Thus, for a > the growth of structures speeds up in comparison 
with the Newtonian case. For the a < case the effect is to diminish the formation of structures. For r > A the 
dynamics is essentially Newtonian. However, this preliminar analysis we have done is insufficient to give us a clear 
conclusions on the role played by SF in the large-scale structure formation process. We will need to do a systematic 
study of the evolution of the two-point correlation function which is a mesure of galaxy clustering. We also will need 
to compute the mass power spectrum and velocity dispersions of the halos. Therefore, we will be able make sistematic 
comparisons with observations. This work is in process and will be published soon. 
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